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Continuous wave gravitational radiation from isolated rotating neutron stars is 
discussed. The general waveform and orders of magnitude for the amplitude are 
presented for various known pulsars. The specific case of gravitational radiation 
resulting from the distortion induced by the stellar magnetic field is presented. 
Finally some preliminary results about the signal from the whole population of 
neutron stars in the Galaxy are discussed. 



1 Introduction 

As compact objects with significant internal velocities, neutron stars consti- 
tute a priori valuable gravitational wave sources. Besides. catastrophic events 
such as their coalescence El or their gravitational collapse □, these stars can be 
interesting sources of continuous wave gravitational radiation, provided they 
deviate from axisymmetry. Whereas the accompanying lecture H focuses on 
mechanisms, such as spontaneous symmetry breaking, which imply accretion 
from a companion, the present lecture is devoted to radiation from isolated 
rotating neutron stars. 

2 General considerations 

2.1 A general formula for gravitational emission by a rotating star 

When dealing with neutron stars, the classical weak- field quadrupole formula 
for gravitational wave generation is not valid. For highly relativistic bodies, 
Ipsera has shown that the leading term in the gravitational radiation held hij 
is given by a formula which is structurally identical to the quadrupole formula 
for weak-field sources u, the Newtonian quadrupole being simply replaced by 
Thome's quadrupole moment U 2^. The non-axisymmetric deformation of 
neutron stars being very tiny, the total Thome's quadrupole can be linearly 
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decomposed into the sum of two pieces: = T, t 
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where X'' ot is the 



quadrupole moment due to rotation and 2^ lst is the quadrupole moment due 
to the process that distorts the star, for example an internal magnetic fieldQ, 
anisotropic stresses from the nuclear interactions or irregularities in the solid 
crustQ. Let us make the assumption that the distorting process has a privileged 
direction, i.e. that two of the three eigenvalues of lfl st are equal. Let then a 
be the angle between the rotation axis and the principal axis of ifj st which 
corresponds to the non-degenerate eigenvalue. The two modes h + and h x of 
the gravitational radiation field in a transverse traceless gauge derived from 
the Thorne-Ipser quadrupole formula arel3 
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where i is the inclination angle of the "line of sight" with respect to the rotation 
axis and 

16tt 2 G Ie 

h ° = —^J*- r > ( 3 ) 

where r is the distance of the star, P = 2tt/H is the rotation period of the star, 
/ its moment of inertia with respect of the rotation axis and e := —3/2X^/1 
is the ellipticity resulting from the distortion process. 

It may be noticed that Eqs. (0)-(||) are structurally equivalent to Eq. (1) 
of Zimmermann & Szedenitsu, although this latter work is based on a different 
physical hypothesis (Newtonian precessing solid star). 

From formula? (|l|)-(||), it appears clearly that there is no gravitational emis- 
sion if the distortion axis is aligned with the rotation axis (a = or 7r) . If both 
axes are perpendicular (a — 7r/2), the gravitational emission is monochromatic 
at twice the rotation frequency. In the general case (0 < \a\ < tt/2), it con- 
tains two frequencies: f2 and 20. For small values of a the emission at 57 is 
dominant. Numerically, Eq. (0) results in the amplitude 
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Note that / = 10 38 kgm 2 is a representative value for the moment of inertia of 
a 1.4M neutron star. In the following, / is systematically set to this value. 

The values of ho resulting from Eq. (Ml) are given in Table for two r-¥Pung 
rapidly rotating pulsars (Crab and Vela), the nearby pulsar GemingalMJ and 
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Table 1: Gravitational wave amplitude ho on Earth as a function of the ellipticity e for five 

pulsars. 



name 


rotation period 


distance 


GW amplitude 




P [ms] 


r [kpc] 


h 


Crab 


33 


2 


1.9 x 10- 2Y (e/10- b ) 


Vela 


89 


0.5 


1.1 x 10- 27 (e/10" 6 ) 


Geminga 


237 


0.16 


4.7 x 10- 28 (e/10~ 6 ) 


PSR B1957+20 


1.61 


1.5 


1.1 x 10- 24 (e/10- 6 ) 


PSR J0437-4715 


5.76 


0.14 


9.1 x 10- 25 (e/10- 6 ) 



two millisecond pulsars: the secondBiastest one, PSR 1957+20, and the nearby 
millisecond pulsar PSR J0437-4715EI. At first glance, millisecond pulsars seem 
to be much more favorable candidates than the Crab or Vela. However, in 
Table [l], e is in units of 10~ 6 and the very low value of the period derivative P 
of millisecond pulsars implies that their ellipticity is at most 2 x 



shall see in S 2.3 



2.2 Detectability by VIRGO 

Let us give a crude estimate of the minimum amplitude h detectable by the 
VIRGO interferometer. Whereas the expected amplitude is very weak, as 
compared with other astrophysical processes such coalescences or gravitational 
collapses, one can take advantage of the permanent character of the signal to 
increase the signal-to-noise ratio by increasing the observing time. Indeed for 
an integration time T, the signal-to-noise ratio reads 

| = J^L^t , (5) 
N h (f) 

where h(f) is VIRGO sensitivity (square root of the noise spectral density) at 
the frequency /. The miminum values of ho leading to S/N = 1 when T = 3 yr 
are given m Table g In view of these values, we shall take as a basis for our 
discussion that in order to be detectable by VIRGO, a rotating neutron star 
must produce ho > 10 -26 . 



'the "historical" millisecond pulsar PSR 1937+21, which is the fastest one, is not consid- 
ered for it is more than twice farther away. 
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Table 2: Minimal amplitude ho detectable (S/N = 1) by VIRGO within 3 mars of integra- 
tion. The values of h(f) have been taken from Giazzoto et alEj 



frequency 
/ [Hz] 


sensitivity 
h [Hz" 1 / 2 ] 


detectable amplitude 
min ho 


10 


10" 2i 


io-' Zb 


30 


1Q -22 


10 -26 


100 


3 x 10" 23 


3 x 10- 27 


1000 


3 x 10" 23 


3 x 10- 27 



Table 3: Absolute upper bounds on the ellipticity and the GW amplitude derived from the 
measured spin-down rate of pulsars. 





GW frequencies 


max. 


max. 


name 






ellipticity 


GW amplitude 




/ [Hz] 


2/ [Hz] 


^max 


^O.max 


Vela 


11 


22 


1.8 x 10" 3 


1.9 x 10" 24 


Crab 


30 


60 


7.5 x 10- 4 


1.4 x 10- 24 


Geminga 


4.2 


8.4 


2.3 x 10- 3 


1.1 x 10- 24 


PSR B1509-68 


6.6 


13.2 


1.4 x 10- 2 


5.8 x 10- 25 


PSR B1706-44 


10 


20 


1.9 x 10" 3 


4.2 x 10- 25 


PSR B1957+20 


621 


1242 


1.6 x 10" y 


1.7 x 10" 27 


PSR J0437-4715 


174 


348 


2.9 x 10" 8 


2.6 x 10- 26 



2.3 Upper bounds on gravitational radiation from pulsars 

An absolute upper bound on ho can be derived by assuming that the observed 
slowing down of the pulsar (the so-called P) is entirely due to the energy car- 
ried away by gravitational radiation. Let us stress that this is not a realistic 
assumption since most of the P is thought to result instead from losses via 
electromagnetic radiation and/or magnetospheric acceleration of charged par- 
ticles — at least for Crab-like pulsars. However, this provides an upper bound 
on the ellipticity e and the gravitational wave amplitude h$. The resulting 
values are given in Table ^. The five first entries in this Table correspond to 
the fpifi-highest values of ho max among the 706 pulsars of the catalog by Taylor 
et alfcM 

New et alta have recently pointed out that if the mean ellipticity of pulsars 
is taken to be of the order of the e max of millisecond pulsars, i.e. e ~ 10~ 9 (cf. 
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Table 4: Minimum values of the ellipticity required to produce ho = 10 on Earth. 



name 


min e 


Crab 
Vela 
Geminga 
PSR B1957+20 
PSR J0437-4715 


5.3 x 10^ = 7xl0~ 3 e max 
9.1 x 1(T 6 =5xl(T 3 e max 
2.1xl(T 5 = 9xl0~ 3 e max 
9.1 x 1(T 9 > e max 
l.lxlCT 8 = 0.4e max 



Table then the Crab pulsar reveals to be a much worse candidate than PSR 
B1957+20, as it can be seen by setting e = 1CT 9 in Table |: /ig lab ~ 2 x KT 30 
versus hT 7+w - lO" 27 . 



2-4 Ellipticity required for a detectable amplitude 

Given the threshold ho = 10 -26 for detectability by the VIRGO interferometer 
(§ 2.2), one can consider the corresponding value of e resulting from Eq. (^) 



and compare it with the maximum value given by the pulsar slowing down 
(Table ||) . The results are presented in Table ||. From that it can be concluded 
that an ellipticity as small as 10~ 8 leads to a detectable amplitude for the 
nearby millisecond pulsar PSR J0437-4715, whereas the Crab or Vela pulsar 
should have an ellipticity of the order 10 -5 , which is about one perc ent of the 
maximum allowable ellipticity as given by the spin-down rate (§ |2.3|) . 



3 The specific case of magnetic field induced distortion 

In this section, we consider the specific example when the distortion results 
from the neutron star's magnetic field. In this case the ellipticity is express- 
ible as e = /3A4 2 /A4q, where M. is the magnetic dipole moment, Aio = 
2.6 x 10 32 Am 2 and (3 a dimensionless coefficient which measures the efficiency 
of the magnetic structure in distorting the star. For an incompressible fluid 
Newtonian body endoved with a uniform magnetic field E£l, (3 — 1/5. For a 
given pulsar, M. can be inferred from the value of PP. 

We have developed a numerical code to compute the deformation of mag- 
netized neutron stars within general relativityQETJ. The solutions obtained are 
fully relativistic and self-consistent, all the effects of the electromagnetic field 
on the star's equilibrium (Lorentz force, spacetime curvature generated by the 
electromagnetic stress-energy) being taken into account. The magnetic field is 
axisymmetric and poloidal. 
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Magnetic field lines 



Magnetic field lines 




Figure 1: Magnetic field lines generated by a current distribution localized in the crust of 
the star (left) or exterior to a type I superconducting core (right). The thick line denotes 
the star's surface and the dashed line the internal limit of the electric current distribution 
(left) and the external limit of the superconducting region (right). The distortion factor 
corresponding to these configurations is j3 = 8.84 (left) and f) = 157 (right). 



The reference (non-magnetized) configuration is taken to be a 1.4M 
static neutron, star built with the equation of state UV14 +TNI of Wiringa, Fiks 
& FabrociniEa. Various magnetic field configurations have been considered; the 
most representative of them are presented hereafter. 

Let us first consider the case of a perfectly conducting interior (normal 
matter, non-superconducting). The simplest magnetic configuration compati- 
ble EJ with the MHD equilibrium of the star results in electric currents in the 
whole star with a maximum value at half the stellar radius in the equatorial 
plane. The computed distortion factor is — 1.01, which is above the 1/5 
value of the uniform magnetic field/incompressible fluid Newtonian model 
but still very low. Another situation corresponds to electric currents localized 
in the neutron star crust only. Figure presents one such configuration: the 
electric current is limited to the zone r > 7% = 0.9 r eq . The resulting distortion 
factor is (3 = 8.84. 

In the case of a superconducting interior, of type I (which means that 
all magnetic field has been expulsed from the superconducting region), the 
distortion factor somewhat increases. In the configuration depicted in Fig. 0, 
the neutron star interior is superconducting up to = 0.9 r eq . For r > r*, the 
matter is assumed to be a perfect conductor carrying some electric current. 
The resulting distortion factor is (3 = 157. For r* ~ 0.95 r cq , (3 is even higher: 
(3 = 517. 

The above values of (3, of the order 10 2 — 10 3 , though much higher than 
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in the simple normal case, are still too low to lead to an amplitude detectable 
by the first generation of interferometric detectors in the case of the Crab 
or Vela pulsar, which would require (3 > 10 4 . It is clear that the more 
disordered the magnetic field the higher /3, the extreme situation being reached 
by a stochastic magnetic field: the total magnetic dipole moment M. almost 
vanishes, in agreement with the observed small value of P, whereas the mean 
value of P 2 _throughout the star is huge. Note that, according to Thompson 
& Duncan E£l, turbulent dynamo amplification driven by convection in the 
newly-born neutron star may generate small scale magnetic fields as strong 
as 3 x 10 11 T with low values of -Bdipoic outside the star and hence a large (3. 
In order to mimic such a stochastic magnetic field, we have considered the case 
of counter-rotating electric currents. The resulting distortion factor can be as 
high as (3 = 5.7 x 10 3 . 

If the neutron star interior forms a type II superconductor, the magnetic 
field inside the star is organized in an array of quantized magnetic flux tubes, 
each tube containing a magnetic field B c ~ 10 11 T. As discussed by Ruderman 
E3, the crustal stresses induced by the pinning of the magnetic flux tubes is of 
the order B c B/2/iq, where B is the mean value of the magnetic field in the crust 
(B ~ 10 8 T for typical pulsars). This means that the crust is submitted to 
stresses ~ 10 3 higher than in the uniformly distributed magnetic field (compare 
B c B/2fio with _B 2 /2^o). The magnetic distortion factor (3 should increase in 
the same proportion. We have not done any numerical computation to confirm 
this but plan to study type II superconducting interiors in a future work. 



4 Gravitational radiation background from the whole population 
of neutron stars in the Galaxy 

4-1 The squared signal 

Let N be the total number of neutron stars in our Galaxy. The response of an 
interferometric detector to the gravitational wave field {h\,h\) [cf- Eqs. (@)- 
(0)] emitted by the i th neutron star is 

h i (t) = Fi(t)h\(t)+F* x (t)h i x (t) , (6) 

where F\ (t) and (t) are beam-pattern factors which depend on the direction 
of the star with respect to the detector arms. They vary with time because of 
the Earth rotation and revolution around the Sun. This results in an amplitude 
modulationEiQ as well as a Doppler shift E2E 3 ] of the signal. The time-average 



of the total (i.e. the sum on all the galactic neutron stars) is zero but not the 
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Figure 2: Total squared signal (divided by N) from a distribution of N neutron stars 
concentrated in the galactic disk with a height scale of 0.5 kpc. The average rotation period 
is assumed to be P = 5 ms and the mean ellipticity t = 10 — 8 . 



time-average of the squared total signal: 

n 2 



M - i / 



io+T 
to 



N 



dt . (7) 



The key-point is that if the galactic neutron star distribution is not isotropic, 
(h 2 ) will exhibit temporal variation, with a period of one sidereal day. The pre- 
cise shape of the signal depends on the neutron star distribution. An example 
lj corresponding to a-disk distribution is presented in Fig. ^. 

It can be shownEJ that the signal-to-noise ratio for detecting this signal is 

S _ {h 2 } ( T V /2 ^ (8) 



N 2h 2 \ 2Aiy 

where Ai> is the frequency bandwith and T is the observation time. It appears 
that such a signal might be detected if the number of neutron stars falling in 
the frequency range of the VIRGO detector is larger than ~ 10 6 - 10 7 . 

4-2 Number of rapidly rotating neutron stars in the Galaxy 

From the star formation rate and the outcome of supernova explosiQiis p _the 
total number of neutron stars (NS) in the Galaxy has been estimatedE3~E3 to 
be about 10 9 . The number of observecUNS-is much lower than this number: 



700 NS are observed as radio pulsars EJ"li-3, ~ 150 as X-ray binaries, among 
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which ~ 30 are_K-ray pulsars ElrEa, and a few as isolated NS, through their 
X-ray emissionLJ. 

For our purpose, the relevant number is given by the fraction of these 
~ 10 9 NS which rotates sufficiently rapidly to emit gravitational waves in 
the frequency bandwith of VIRGO-like detectors. Putting the low frequency 
threshold of VIRGO to / min = 5 Hz (optimistic value !), and using the fact that 
the highest gravitational frequency of a star which rotates at the frequency / 



is 2f (cf. § 2.1), this means the rotation period of a detectable NS must be 
lower than P max = 0.4 s. The NS that satisfy to this criterion can be divided 
in three classes: (CI) young pulsars which are still rapidly rotating (e.g. Crab 
or Vela pulsars); (C2) millisecond pulsars, which are thought to have been 
spun up by accretion when being a member member of a close binary system 
(during this phase the system may appear as a low- mass X-ray binary); (C3) 
NS with P < 0.4 s which do not exhibit the pulsar phenomenon. 

The number of millisecond pulsars in the Galaxy is estimated E3 to be of 
the order N% ~ 10 5 . The number of observed millisecond pulsars (P < 10 s) is 
about 50 and is continuously increasing. 

The number of young (non-recycled) rapidly rotating NS is more difficult 
to evaluate. An estimate can be obtained from the fact that the observed non- 
millisecond pulsars with P < 0.4 s represent 28 % of the number of catalogued 
pulsars and that the total number of active pulsars in the Galaxy is around 
5 x 10 5 . The number of rapidly rotating non-recycled pulsars obtained in this 
way is Ni ~ 1.4 x 10 5 . 

Adding N± and N2 gives a number of ~ 2 x 10 5 NS belonging to the 
populations (CI) and (C2) defined above. This number can be considered as 
a lower bound for the total number of NS with P < 0.4 s. The final figure 
depends on the number of members of the population (C3) . This latter number 
is (almost by definition !) unknown. All that can be said is that it is lower 
than the total number of NS in the Galaxy (~ 10 9 ). 
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